#include "QR.h"

using namespace mylib;
int main(int argc, char const *argv[])
{
  double ma[] = {17, 24, 1, 8, 15,  //
                 23, 5, 7, 14, 16,  //
                 4, 6, 13, 20, 22,  //
                 10, 12, 19, 21, 3, //
                 11, 18, 25, 2, 9};
  eigen<double> e(5, ma);
  e.HessenbergA.show(5, 5);
  e.Householder();
  e.HessenbergA.show(5, 5);
  for (int i = 0; i < 50; i++)
  {
    e.shift(-e.HessenbergA[e.n-1][e.n-1]);
    e.Givens();
    e.shift(e.HessenbergA[e.n-1][e.n-1]);
  }
  e.HessenbergA.show(5, 5);
  return 0;
}
